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Abstract. For the nonthermal radio emission of the Galactic Center Arc in situ electron acceleration is imperative. The observed 
radio spectrum can be modeled by a transport equation for the relativistic electrons which includes particle acceleration by 
electric fields, momentum diffusion via scattering by magnetohydrodynamical turbulence and energy losses by synchrotron 
radiation. The accelerating electric fields can be regarded as a natural consequence of multiple reconnection events, caused by 
the interaction between a molecular cloud and the Arc region. The radio spectrum and even the recently detected 150 GHz 
emission, explicitely originating from the interaction regions of a molecular cloud with the magnetized Arc, can be explained 
in terms of quasi-monoenergetically distributed relativistic electrons with a typical energy of about 10 GeV accelerated in 


stochastically distributed magnetic reconnection zones. 
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1. Introduction 


The Galactic Center (GC) Arc is a unique nonthermal radio 
continuum structure in the Galaxy. It is located at a projected 
distance of about 30 pc from the Sgr A complex and consists 
of a network of magnetic filaments filled with relativistic elec- 
trons running almost perfectly perpendicular to the Galactic 
plane (e.g. Yusef-Zadeh et al. 1984). The nonthermal nature 
of these filaments has been proven beyond any doubt by radio 
emission which exhibits a degree of linear polarization which 
is close to the intrinsic value at high frequencies of 60% (Reich 
et al. 1988; Lesch & Reich 1992). We emphasize the nonther- 
mal character of the Arc since its radio spectrum is a very 
remarkable one. A decomposition of its radio spectrum be- 
tween 843 MHz and 43 GHz by (Reich et al. 1988) shows a 
spectral index a = +0.3 (we use the conventional relation that 
the observed flux S, scales as v”). This finding of an increas- 
ing radio spectral index has been confirmed by high resolution 
VLA-observations (Anantharamaiah et al. 1991). 

Such an inverted radiation spectrum with an index of +0.3 
is expected from a quasi-monoenergetic electron distribution 
or an energy distribution with a well-defined low energy cut- 
off, respectively (Lesch et al. 1988). Observations at high fre- 
quencies (between 32 GHz, Lesch & Reich 1992, and 43 GHz, 
Sofue et al. 1987b) seem to indicate a spectral turn over, 
i.e. a fading of the Arc towards higher frequencies. However, 
Reich et al. (2000) detected an enhanced emission at 150 GHz 
slightly offset relative to the most intense vertical nonther- 
mal filaments seen at lower frequencies. This emission origi- 
nates from the apparent interacting areas of dense molecular 
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material with the Arc. 150 GHz emission was also detected 
south of a molecular cloud, where the vertical filaments of the 
Arc cross a weak filamentary structure. The spectrum of the 
emission is inverted relative to 43 GHz and compatible with an 
origin from quasi-monoenergetic electrons or an electron dis- 
tribution with a low energy cutoff, but not with optically thin 
emission from cold dust. Reich et al. (2000) conclude “that the 
coincidence of enhanced emission with regions of interacting 
molecular gas strongly suggests that high-energy electrons are 
accelerated in those places where the magnetic field is com- 
pressed”. They calculated a Lorentz factor of the electrons to 
be y ~ 2 x 10+ emitting synchrotron radiation at 150 GHz 
within a magnetic field of 1 mG. 

To summarize, in the Galactic center several areas appear 
to be filled with relativistic electrons whose energy must be 
distributed quasi-monoenergetically around a few GeV. We 
note that even the very center of our Galaxy Sgr A* ex- 
hibits a radio to infrared spectrum which is in surprisingly 
good agreement with optically thin synchrotron emission of a 
quasi-monoenergetic electron distribution with a typical energy 
of 80 MeV, i.e. Lorentz factor of 160 (Duschl & Lesch 1994). 

In a former paper about particle acceleration in the Arc 
(Lesch & Reich 1992) we considered the magnetohydrody- 
namical interaction of a molecular cloud with the Arc fila- 
ments. Especially we investigated the role of the moving gas 
cloud as a trigger mechanism for magnetic field amplification 
accompanied by dissipation of the magnetic energy via mag- 
netic reconnection. 

Magnetic reconnection takes place when magnetic field 
lines with antiparallel directions encounter. Such a situation 
corresponds to the formation of an electric current sheet. 
Typically astrophysical plasmas are ideal electric conductors, 
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i.e. their electrical conductivity is very high. In such media the 
magnetic field is frozen into the motion of the conducting fluid. 
Any plasma velocity distortion onto the magnetic field lines is 
automatically related to an electric field which is perpendicular 
to the plasma velocity and the magnetic field. This electric field 
is described by the ideal Ohm’s law 


1 
E+-v0xB=0. (1) 
C 


In a plasma where the magnetic field is strongly distorted by 
shear flows, radial explosive flows or stochastic motions, it 
is unavoidable that field lines with antiparallel directions en- 
counter. Consequently, the ideal conducting plasma switches 
into a nonideal medium with a finite, localized electrical con- 
ductivity, i.e. the ideal form of Ohm’s law is violated. The 
strong spatial gradients of the magnetic field in such interaction 
zones represent an energetic “crisis” which is relaxed by partial 
dissipation of the magnetic energy via the formation of current 
sheets in the resistive medium (e.g. Priest & Forbes 2000, and 
references therein). The nonideal Ohm’s law is given by 


1 
E+—vxB=R#0, (2) 
C 


where R is some yet unspecified nonideal term, i.e. the plasma 
resistance. 


In other words, magnetic reconnection corresponds to a 
magnetic field aligned electric field which is associated with 
a generalized electric potential (Schindler et al. 1988, 1991) 


v=- fE ds=- [Ras 


where the integral is evaluated along the magnetic field lines 
that penetrate the reconnection region. Of course, such a po- 
tential drop offers the possibility to accelerate particles very ef- 
ficiently (Schindler et al. 1991; Blackman 1996; Lesch & Birk 
1997; Lesch & Birk 1998; Schopper et al. 1998; Litvinenko 
1999; Nodes et al. 2003). 


Instead of investigating the acceleration in one reconnec- 
tion zone, we consider in this contribution the acceleration of 
relativistic electrons in numerous reconnection regions, i.e. cur- 
rent carrying filaments driven by the interaction of a molecular 
cloud with the magnetic field in the Galactic Center Arc (Lesch 
& Reich 1992; Serabyn & Morris 1994). This scenario arises 
quite naturally, since the necessary energy source is represented 
by moving molecular gas which encounters the poloidal mag- 
netic fields in the Arc which is proven to be present by the ob- 
served very high polarization of the radio emission up to 60%. 


(3) 


Since we have many acceleration regions we apply a dy- 
namical description of the energy distribution function of rel- 
ativistic electrons based on systematic momentum gains by 
reconnection, losses by synchrotron losses and momentum 
diffusion due to scattering on magnetohydrodynamical turbu- 
lence, i.e. Alfvén waves and magnetosonic waves. 
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2. The temporal evolution of the energy 
distribution function of relativistic particles 


The dynamical evolution of a population of high-energy par- 
ticles that do not interact with each other is described by the 
Liouville equation 

OF dr d 

+ VrF + = VpF =0 
where F(r, p, t) is the one-body distribution function. Due to 
the complexity of the non-linear dynamics in 6+1 dimensions 
this equation is only of limited use. Fortunately, for our pur- 
poses we do not need to know the full information of the 
particle dynamics. Rather, we will make use of the follow- 
ing assumptions. First, we will assume pitch angle isotropy. 
This should be granted by efficient particle scattering in mo- 
mentum space in reconnection regions caused by Alfvénic and 
magnetosonic wave fields with relatively high-energy densities 
(see Schlickeiser 1986). Second, we are interested in the global 
energy spectrum of the radiating electrons and thus, can dis- 
pense with the detailed spatial dependence of the distribution 
function. 

Consequently, we deal with the isotropic distribution 
f(po = f drF (r, p,t)/V where V denotes the emission 
volume. The total particle number is calculated as N(t) = 
An f fp°dp. The continuous momentum gains and losses Ap 
of the particles can be described by the momentum operator 


(4) 


Lp g E [ Peain F pios)p’] > (5) 
where the systematic temporal momentum gains and losses are 
denoted by Pgain and Pioss, respectively. 

Additionally, we will allow for stochastic diffusion in mo- 
mentum space by scattering of magnetohydrodynamical fluc- 
tuations (in the sense of a Fokker-Planck term) denoted by D,. 
In principle, catastrophic losses that may result to a further sink 
in the balance equation for f can be modeled by a character- 
istic time T, i.e. 0f/Ot ~ —f/t. Finally, the injected particle 
population is represented by some explicit source term S(p). 
Accordingly, the temporal evolution of the particle distribution 
is governed by the equation (Schlickeiser 1984, 1986) 


Of(p,t) | ale 2 Of (p, t) fp D 
ôt p2op\ ? Op T 


J+ Lor (p, t)+ =S (pX(6) 
which holds, if the distribution function f(r, p, t) is separable in 
space and momentum, and if the momentum operator Lp does 
not depend on the spatial coordinates. In the present application 
based on the reconnection scenario (see also Birk et al. 2001) 
the momentum gain is given by Pgain = 4 (ymv) = qE,=f= 
const., where y denotes the Lorentz factor, E, is the acceler- 
ating electric field component (locally directed parallel to the 
magnetic field) and q is the electrical charge. 

We emphasize that the acceleration of high-energy parti- 
cles is regarded as a consequence of the conversion of mag- 
netic field energy into particle energy by magnetic reconnec- 
tion. Since synchrotron losses are responsible for the observed 
non-thermal radio emission from the Arc we take Doss = —kp* 
with k = to7Up Ime, with the Thomson cross section or, 
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the electron mass me, the speed of light c and the magnetic 
field energy density Ug. For our applications inverse Compton 
scattering, which has the same y-dependence as synchrotron 
losses, as well as bremsstrahlung are negligible. In case of par- 
ticle scattering in momentum space by Alfvénic waves the mo- 
mentum diffusion coefficient is given (Schlickeiser 1984) by 
Dy = Dp? with D = v,/9K(p), where vag is the Alfvén veloc- 
ity and K(p) = AcB/3 is the spatial diffusion coefficient, which 
has been expressed in terms of the particle’s mean free path 2 
(Schlickeiser 1986). For our purpose we have set 8 = v/c = 1. 
Since we have to model the formation of monoenergetic distri- 
bution functions, the effect of catastrophic losses must be neg- 
ligible, because otherwise the particles would leave the accel- 
eration region before a pile up in energy has been established 
(Schlickeiser 1986). 


3. Numerical results 


In order to calculate the temporal evolution of f(p, t) by numer- 
ical integration of Eq. (6) we first have to specify the physical 
parameters ¢, k, D, and S (p), which are related to the physical 
conditions of the Galactic Center Arc. First of all, we discuss 
the source term S(p). Since the energy distribution of the in- 
jected particles, which are generated in the Galactic center, is 
a monoenergetic one (Lesch & Reich 1992), the use of a 6- 
like source term, represented, e.g., by a cosh“*(p)-function in 
our calculations, is a reasonable assumption. As a typical injec- 
tion energy we use 75 MeV, the value close to the one calcu- 
lated by (Duschl & Lesch 1994) for the central object Sgr A*, 
i.e. Yinj = 150. 

The magnetic field strength in the arcs filaments is of the 
order of a few mG, or higher (Yusef-Zadeh & Morris 1987a). 
To be on the save side we use the 150 GHz detection of (Reich 
et al. 2000) close to the Arc as the maximum synchrotron fre- 
quency emitted. Thus, we have to deal with a Lorentz factor of 
about y ~ 2x 10*, necessary for 150 GHz synchrotron emission 
in a magnetic field of 107° G. 

The maximum electric field strength is given by E ~ w/cB, 
where w denotes the bulk velocity of the thermal plasma that 
acts as an MHD generator. The detected cloud motion is in the 
range of 15—45 kms‘! (Tsuboi et al. 1997). Again to be on the 
save side we will use v = 5 kms~! which gives a maximum 
electric field strength of 1078 statvoltcm™!. If we assume, that 
electrons can be accelerated over the whole extension of the 
Arc (about 40 pc) without energy losses (ideal linear accelera- 
tor), we need a minimum strength of the average electric field 
of Es ~ 3 x 107! statvoltcm™! to reach the observed particle 
energies, which is a tiny fraction of the available electric field. 

This oversimplification is helpful to work out the analysis 
with a reference value of ¢ and is not meant as a physical as- 
sumption we apply (see discussion below). 

If we restrict our calculations to electrons we obtain an es- 
timate for the systematic acceleration term 
Pgain = CEs = £ = 1.5 X 10” dyne. (7) 
Due to the high polarization of the synchrotron radiation, up 
to 60% (Lesch & Reich 1992) we consider only synchrotron 
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emission to be responsible for the energy losses of the par- 
ticles in the Arc. In this case we can find a certain value 
po = 3ém2c?/4c07Ug)!/* where momentum gains and losses 
keep the balance. For momenta p > pọ momentum loss domi- 
nates momentum gain. For momenta p < po the opposite holds. 
Thus, we expect a pile up for the distribution function f(p, t) 
at po, which is consistent with a monoenergetic distribution 
of the electrons (see also Lesch & Reich 1992). As discussed 
above, the pile up in the Arc should occur at momenta with 
Lorentz factor y ~ 2 x 10+. If one assumes, for simplic- 
ity, a constant electric field of Es ~ 3 x 107! statvolt cm! 
over the complete expansion of the arc, a magnetic field 
of 3 x 107? Gauss is necessary for a pile up at po = ymec = 
5.4 x 107! gcms™! according to a systematic loss coefficient 
of 

pote ease g! cm™!. 


3 me 


(8) 


However, the particles, in fact, will be accelerated locally in nu- 
merous reconnection regions with E£ + 0 and will be scattered 
by magnetohydrodynamic fluctuations like Alfvénic and/or 
magnetosonic waves. Consequently, the electric field strength 
in the reconnection regions must be ultimately stronger 
than 107! statvoltcm™!. Indeed, the electric field can be es- 
timated by the ratio of the number density of relativistic elec- 
trons n, and the thermal electron density ne (Papadopoulos 
1977) 
nr E. 
T ~ exp|->=| (9) 
Ec = MeVtheVcou /e denotes the critical electric field, whereas Vihe 
is the thermal velocity of the electrons and veon is the colli- 
sion frequency. As calculated by Lesch & Reich (1992) the 
thermal electron number density ne in the Arc filaments is 
about 10 cm~? and the density of relativistic electrons is of the 
order of 10-7 cm~?. If we use an electron temperature of 10* K 
we obtain from Eq. (9) E = 3 x 107!° statvoltem™!, which can 
be easily achieved by the available electric field, induced by the 
w x B -motion of a molecular clouds calculated above. 

To obtain the diffusion coefficient D we first have to esti- 
mate the particle’s mean free path 2. As shown in (Schlickeiser 
2002) A can be determined by 


(10) 


where rg is the gyro radius of the protons, Bp /6B gives the ratio 
between macroscopic and fluctuated magnetic field strength, 
Lo is the expansion of the considered object and q is turbulence 
spectral index. In the vicinity of a strong magnetic distortion, 
like the Arc filaments, Burgers turbulence is expected to be ex- 
cited (Chambers et al. 1988). Thus, we choose a turbulence 
spectral index of q = 2. With a particle density of n ~ 10 cm™ 
and a magnetic field strength of 3x 107° Gauss the Alfvénic ve- 
locity is about 700 km s~!. Hence we obtain for the diffusion 
coefficient 

vy 1 (5B) v -18 ei 

p= di =3(F)  -saxio Ss (11) 
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if we assume, that Bo/ôB is of the order of 10 which im- 
plies a reasonable level of well-developed turbulence. Much 
higher diffusion coefficients would inhibit pile-up distributions. 
However, a quasi-linear theory applied for Eq. (9) requires 
OB < Bo. 

The synchrotron emission spectrum associated with the 


particle spectrum can be calculated from (Rybicki & Lightman 
1979) 


fd f P,(p)N(p) dp, (12) 
where 

2 3 Dy); fore) 
Pp) = a i „KOR, (13) 


is the power per unit frequency v emitted by each electron. 
K ; (€) denotes the modified Bessel function of the second kind 


and ve = 3y3we/4r is the critical frequency (we is the electron 
gyro frequency). 

Having clarified the initial conditions and chosen physical 
parameters we now come to the numerical results. Figure 1 
shows the temporal evolution of the particle density N(y) and 
the corresponding synchrotron emission spectrum, calculated 
from Eq. (12). In the chosen normalization the time is mea- 
sured by fọ = 10° s. As long as the particle acceleration ex- 
ceeds the synchrotron losses, the quasi-monoenergetic distri- 
bution of the electrons is shifted toward higher momenta (a, c). 
Due to momentum diffusion the maximum of N(y) decreases 
with time (c). As soon as the synchrotron losses are compara- 
ble to the energy gains, the maximum of N(y) increases again 
(Fig. 2a). At a Lorentz factor of about y = 2 x 104 the momen- 
tum gains and losses are balanced and the distribution of the 
electrons that are assumed to be contained in the Arc during the 
entire radiation processes develops towards a shifted pile up. 
Consequently, a monoenergetic distribution function evolves 
again (Fig. 2c). The cut-off frequency remains almost constant 
during entire pile-up dynamics. The associated radiation spec- 
trum (Fig. 2c) corresponds to the one observed by several au- 
thors, including the 150 GHz emission detected by (Reich et al. 
2000). We note that for a finite loss time of r < 10!” s the pile- 
up is prevented by the particle losses. 


4. Conclusions 


In the filaments of the Galactic Center Arc electrons are accel- 
erated to considerably high energies of about 10 GeV. This is 
obvious from the detected polarized radio emission exhibiting 
a very high degree of polarization (60%) and a rising spectrum 
with a spectral index of +0.3 up to 150 GHz. Such a spectral 
behavior is due to an energy distribution function of the radi- 
ating relativistic electrons which either has a low-energy cutoff 
or which is quasi mononenergetic. No obvious energy sources 
for particle acceleration are present in the neighborhood of the 
filaments. The very center of the galaxy, SgrA* is a source for 
monoenergetic relativistic particles but on significant lower en- 
ergies of about 50 MeV (Duschl & Lesch 1994). However, 
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a) Temporal evolution of N(y). The solid, dotted 
and dashed lines show the injected particle spec- 
trum and the evolution after t = 1 tọ and t = 2 tọ, 
respectively. 
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b) The corresponding synchrotron spectrum 7,. 
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c) Temporal evolution of N(y). The solid and dot- 
ted lines show the evolution after t = 10 tọ and t = 
20 to, respectively. 
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Fig. 1. Temporal evolution of N(y) calculated from Eq. (6). The val- 
ues of the physical parameters, which we have used for our calcu- 
lations are: ¢ = 1.5 x 107” dyne, k = 1.4 x 10° g-'cm™! and 
D = 5.4 x 107!8 s“!. To evaluate the synchrotron emission spectrum 
(see Eq. (12)) we have used a magnetic field strength of B = 3 mG. 
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a) Temporal evolution of N(y). The solid and dot- 
ted lines show the evolution after t = 30 tọ and 
t = 50 fo, respectively. 
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b) The corresponding synchrotron spectrum /,. 
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c) Temporal evolution of N(y). The solid, dotted 
and dashed lines show the evolution after t = 80 fo, 
t = 100% and t = 150 tọ, respectively. 
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Fig. 2. Temporal evolution of N(y) towards a quasi-monoenergetic 
distribution function a), c). The resulting synchrotron emission spec- 
trum d) shows an inverted spectrum with a spectral index of œ = 
1/3 Ud, ~ v”) and a cutoff-frequency of some hundred GHz af- 
ter 1.5 x 10!! s. 
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the magnetic filaments interact with the plasma of a molecu- 
lar cloud which moves with velocities of some 10 kms"! rela- 
tive to the magnetic field. This interaction can be interpreted in 
terms of the induced Lorentz force inducing a convective elec- 
tric field which in the first place is oriented perpendicular to 
the magnetic field and the local plasma velocity of a molecular 
cloud. The motion of the plasma locally distorts the magnetic 
field necessarily in such a way that antiparallel directed field 
lines encounter. Provided that some violation of ideal Ohm’s 
law occurs, e.g. due to microturbulence, magnetic reconnec- 
tion will occur. The magnetic field energy will be partly con- 
verted to particle energization. In the reconnection regions a 
magnetic field-aligned electric field component forms with a 
magnitude of some fraction of the convective electric field. This 
parallel electric field represents a perfect candidate for efficient 
particle acceleration. The required electric field for the elec- 
tron energization are considerably weaker than the magnetic 
field which is in accordance with the observational fact that 
the cloud plasma moves with velocities much lower than the 
speed of light. Since reconnection is a localized phenomenon 
which depends on the local properties of the plasma, we con- 
sider a scenario characterized by multiple reconnection sheets 
in which particles are accelerated and scattered by magnetohy- 
drodynamical fluctuations. Such scattering leads to diffusion in 
energy space. Moreover, we take into account the energy losses 
by the observed synchrotron radiation. By means of a relatively 
simply transport equation for the distribution function of the ra- 
diating relativistic electrons we could show that an injected mo- 
noenergetic distribution function with 75 MeV (coming from 
the Galactic center) evolves into a quasi-monoenergetic distri- 
bution function at 10 GeV. The resulting radiation spectra in- 
crease with frequency up to some hundred GHz according to 
the observations, i.e. S, ~ v!/3 

A major assumption of our study is the complete isotropy 
of the distribution function which is based on calculations 
by (Achatz et al. 1990), who show that the isotropisation 
time fiso = c/va(B/ôB)1/Q, where Q is the gyro frequency of 
the protons. In this application fiso is about 4 x 10° s which 
is much shorter than the studied time scale of the evolution of 
the distribution function. Some non-isotropic particle popula- 
tion would lead to an overlap of a second radiation component 
in a specific energy range that we can not handle in our model. 
However, the present observations seem to give no hint to such 
a component. 
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